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ABSTRACT 



o 



Context. Knowledge of the intrinsic shape of galaxy components is a crucial piece of information to constrain phenomena driving 
their formation and evolution. 

Aims. The structural parameters of a magnitude-limited sample of 148 unbarred SO-Sb galaxies were analyzed to derive the intrinsic 
shape of their bulges. 

Methods. We developed a new method to derive the intrinsic shape of bulges based on the geometrical relationships between the 
apparent and intrinsic shapes of bulges and disks. The equatorial ellipticity and intrinsic flattening of bulges were obtained from the 
length of the apparent major and minor semi-axes of the bulge, twist angle between the apparent major axis of the bulge and the 
galaxy line of nodes, and galaxy inclination. 
■ Results. We found that the intrinsic shape is well constrained for a subsample of 1 15 bulges with favorable viewing angles . A large 

fraction of them is characterized by an elliptical section (B/A < 0.9). This fraction is 33%, 55%, and 43% if using their maximum, 
mean, or median equatorial ellipticity, respectively. Most are flattened along their polar axis (C < (A + B)/2). The distribution of 
triaxiality is strongly bimodal. This bimodality is driven by bulges with Sersic index n > 2, or equivalently, by the bulges of galaxies 
with a bulge-to-total ratio B/T > 0.3. Bulges with n < 2 and with B/T < 0.3 follow a similar distribution, which is different from 
that of bulges with n > 2 and with B/T > 0.3. In particular, bulges with n < 2 and with B/T < 0.3 show a larger fraction of 
oblate axisymmetric (or nearly axisymmetric) bulges, a smaller fraction of triaxial bulges, and fewer prolate axisymmetric (or nearly 
axisymmetric) bulges with respect to bulges with n > 2 and with B/T > 0.3, respectively. 

Conclusions. According to predictions of the numerical simulations of bulge formation, bulges with n < 2, which show a high 
' fraction of oblate axisymmetric (or nearly axisymmetric) shapes and have B/T < 0.3, could be the result of dissipational minor 

mergers. Both major dissipational and dissipationless mergers seem to be required to explain the variety of shapes found for bulges 
Tj" \ with n > 2 and B/T > 0.3. 
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1. Introduction 1.1. Intrinsic shape of elliptical galaxies 

The first attempt t o derive the int rinsic shape of elliptical galax- 

rr U u , r i j j i „ 11 i • 1 • 1 ■ ies was done by Hubble (1926). The distribution of their in- 

The halos of cold dark matter assembled in cosmological sim- . . _ . J ' f \ — r\. , , 

«_] ■ , .• . , , . • • , , rrn — j T 1 Ihnnd tnnsic rlattemngs was obtained from the observed ellipticities 

H , ulations appear to be strongly triaxial (see Allgood et al. 2006, , .,„••,,. 

, and references therein). Their intrinsic shape is characterized by ™der the assumption that elliptical galaxies were oblate el- 

■ • . . , • .. D/ , , , . . , • lipsoids with a random orientation with respect to the line of 

an intermediate-to-long axis ratio B/A and a short-to-long axis . F , _ , ,. ., . 

ratio CI A which can vary as a function of radius. On the con- sl % hL Earl y studies considered elliptical galaxies to be axisym- 

trary, the halo shape inferred from observations of the Milky metric system s Oblateness andDj^toiess were assumed by 

(Oiling & Merrifield 2000; Ibata et al. 2001; Johnston et al] I Sandage et al.| m$ and|BinneyJ (EH!, respectively to repro- 
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Way (. ... „ ..„.....^.„ _ 

innci 7" ~u' e^~A-"^"\ u* " ~ as — j (j uce the distribution of observed ellipticities of the Reference 

2005) and a number of individual nearby galaxies (Memfielc _ , 1 , T , L — : „ . ,, : 1 

7004^ i« „«,rlv airisvmmrfTir Th P «t„Hv nf A.. infrin«V s b fl n P Catalog of Bright Galaxies (de Vaucouleu rs & de Vaucou leurs 



20041) is nearly axisymmetric. The study of the intrinsic shape FTq^ ^ ^ j^p^)^ 



of the luminous galactic components may serve to constrain the 

halo shape, which is related to the final morphology of the galaxy Afterwards, a number of kinematic and photometric find- 

and depen ds on the phenome na driving its formation and evolu- ings led to the suggestion that there are also elliptical galaxies 

tion (e.g.. IHeller et al.l 120 07). The intrinsic shapes of elliptical with a triaxial shape. In fact, the low ratio between rotational 

galaxies and disks were extensively studied in the past, whereas velocity and velo city dispersion (Be rtola & Capacciolil 1975: 

bulges appear to be less studied, even if the y account for abou t rillingwort hlll977h, the twisting in the i sophotes (ICarten Tl978; 

25% of the stellar mass of the local universe (IDriver et a!1l2007l) . Bertola & Gallettal [19791 iGallettal Tl 9 8 ) . and the rotation mea- 
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sured along the minor axis (ISchechter & Gunnlll979h of some 
elliptical galaxies could not be exp lained in terms of axisymmet- 
ric ellipsoids. As a consequence. iBenacchio & Gallettal (l980) 
and [Binnev & de Vaucouleursl d 19811) showed that the distribu- 
tion of observed ellipticities could be satisfactorily accounted for 
also in terms of a distributi on of triaxial e l lipsoid s . Similar con- 
clusions were reached bv [Fasano & Vio l (1 19911) . lLambas et alJ 
dl992l) . iRvdenl d 19921 1 19961) . and lFasanol d 19951) . However, dif- 
ferent galaxy samples and different assumptions on triaxiality 
resulted in different distributions of intrinsic axial ratios. In ad- 
dition, not all the ellip tical galaxi e s hav e the same intrinsic 
shape. In fact, Tremblay & Merritt] d!996l) found that the dis- 
tribution of the observed ellipticities of galaxies brighter than 
Mb — -20 is different than that of the less luminous ones. In 
particular, there is a relative lack of highly-flattened bright ellip- 
ticals. This reflects a difference in the shape of low-luminosity 
and high-luminosity ellipticals: fainter ellipticals are moderately 
flattened and oblate , while brighter el lipticals are rounder and 
triaxial. Recently, iFasano etal found also that even if 

both normal ellipticals and brightest cluster galaxies (BCG) 
are triaxial, BCGs tend to have a more prolate shape, and that 
this tendency to prolateness is mainly driven by the cD galax- 
ies present in their sample of BCGs. These kinds of statisti- 
cal analyse s benefit from lar ge galax y samples, such a s thos e 
studied by iKimm & Yil d2007l) and iPadilla & Strauss! d2008h . 
They analyzed the observed ellipticities of early-type galaxie s 
in Sloan Digital Sky Survey dAdelman-McCarthv et~alll2006b . 
Furthermore, these large datasets allowed them to study the de- 
pendence of the intrinsic shape on other galaxy properties, such 
as the luminosity, color, physical size, and environment. 

Determining the distribution of the intrinsic shape of ellipti- 
cal galaxies is also possible by combining photometr ic and kine- 
matic information dBinnevlll985l: iFranx et al.lll991l) . However, 
the resulting distribution of intrinsic flattenings, equatorial el- 
lipticities, and intrinsic misalignments between the angular mo- 
mentum and the intrinsic short axis can not be derived uniquely. 
Only two observables are indeed available, the distribution of 
observed ellipticities and the distribution of kinematic misalign- 
ments between the photometric minor axis and the kinematic 
rotation axis. Therefore, further assumptions about the intrin- 
sic shape and direction of the angular momentum are needed to 
simplify the problem. In addition, this approach requires a large 
sample of galaxies for which the kinematic misalignment is ac- 
curately measured. But, to date this informat ion is available only 
for a few tens of galaxies dFranx et al.ll 19911) . 

Many individual galaxies have been investigated by de- 
tailed dynamical modeling of t he kinematic s of gas, stars, 
and pl anetary nebu la e (e.g ., Tenies etal] 



1994 IStatler & Frvl 1 1994 iM athieu & Deionghe 



2007k iThomas et al. 



20071: Ide Lorenzi et alJ 120091) . Recently. 
van den Bosch & van de Vera d2009l) have investigated how 




Gerhard et al] 1200 lk iGebhardt et all 12003k ICappellari et alJ 



well the intrinsic shape of elliptical galaxies can be recovered 
by fitting realistic triaxial dynamical models to simulated 
photometric and kinematic observations. The recovery based on 
orbit-based models and state-of-the-art data is degenerate for 
round or not-rotating galaxies. The intrinsic flattening of oblate 
ellipsoids is almost only constrained by photometry. The shape 
of triaxial galaxies is accurately determined when additional 
photometric and kinematic complexity, such as the presence 
of an isophotal twist and a kinematically decoupled core, is 
observed. Finally, the intrinsic shape of individual galaxies can 
be also constrained from the observed ellipticity and isophotal 



twist by assuming the intrins ic density distribution dWilliamsl 
119811: [Chakraborty e t al. 200^). 



1.2. Intrinsic shape of disk galaxies 

Although the disks of lenticular and spiral galaxies are often 
considered to be infinitesimally thin and perfectly circular, their 
intrinsic shape is better approximated by flattened triaxial ellip- 
soids. 

The disk thickness can be directly determined from 
edge-on galaxies. It depends both on the wavelength at 
which disks are observed and on galaxy morphological 
type. In deed, galactic disks become t hicker at longer wave- 
lengths dDalcanton & Bernstein! l2002k iMitronova et al] l2004l) 
and late-type spirals have thinner d isks than early-type spirals 
dBottinelli et allll 1 983k [Guthridl 1 992h . 

Determining the distribution of both the thickness and ellip- 
ticity of disks is possible by a statistical analysis of the distribu- 
tion of apparent axia l ratios of randomly oriented spiral galaxies. 
ISandage et al .1 d 19701) analyzed the spiral galaxies listed in the 
RC1. They concluded that disks are circular with a mean flat- 
tening {CI A) - 0.25. However, the lack of nearly circular spiral 
galaxies {B/A 1) rules out that disks h a ve a perfectly axisym- 
metric sha pe. Indeed. | B inggelil dl980h. IBenacchio & Ga lletta 
( 1980), and Bi nnev & de Vaucouleursl d 19811) showed that disks 
are slightly elliptical with a mean ellipticity (1 - B/A) = 0.1. 
These early findings were based on the analysis of photographic 
plates of a few hundreds of galaxies. Later, they were con- 
firmed by measuring ellipticities of several thousands of ob- 
jects in CCD images and digital scans of plates obtained in 
wide-field surveys. The large number of objects allows the con- 
straining of the distribution of the intrinsic equatorial elliptic- 
ity, which is well fitted by a one-sided Gaussian centered on 
1 - B/A = with a standard deviation ranging from 0.1 to 
0.2 and a mean of 1 dLambas et al] 1 1992k iFasanoet al.ll 1993k 
lAlam & Rvde"iill2002k lRvdenl l2004). Like the flattening, the in- 
trinsic ellipticity depends on the morphological type and wave- 
length. The disks of early-type spirals are more elliptical than 
those of late-type spirals an d their media n ellipticity increases 
with observed wavelength dRvdenl 120061) . Furthermore, lumi- 
nous spiral galaxies tend to have th icker and rounder disk s 
than low-luminosity spiral galaxies (Padilla & Strauss 2008). 
Different mechanisms have been proposed to account for disk 
thicken ing, including the scattering of s tars off giant mo lecular 
clouds dSpitzer & Schwarzschiidlll95 lk IVillumsenlll985h tran- 
sient density waves of the sp iral arms (Barbanis & Woitien l967: 
ICarlberg & Sellwoodlll985[). and minor mergers w ith satellite 
galaxies (e.g JOuinn et alJll993klWalker et allll996l) . 

The study of the intrinsic shape of bulges presents similar- 
ities, advantages, and drawbacks with respect to that of ellip- 
tical galaxies. For bulges, the problem is complicated by the 
presence of other luminous components and requires the iso- 
lation of their light distribution. This can be achieved by per- 
forming a photometric decomposition of the galaxy surface- 
brightness distribution. The galaxy light is usually modeled as 
the sum of the contributions of the different galactic compo- 
nents, i.e., bulg e and disk, and eventually lenses, ba rs, spiral 
arms, and rings dPrieto et al.ll2"00l:lAguerr i et al. 2Q05|). A num- 
ber of two-dimensional parametric decomposition techniques 
have been develop e d in the last several years with this aim 
(e.g.. ISimardlll998k iKhosroshahi et al.ll2000k iPeng et al.l r2 002: 
de Souza et al. 2004; Laurikaine n et al] 12005k iPignatelli et al.l 
2006; iMendez-Abreu et alj|2008l) . On the other hand, the pres- 
ence of the galactic disk allows for the accurate constraining of 
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the inclination of the bulge under the assumption that the two 
components share the same polar axis (i.e., the equatorial plane 
of the disk coincides with that of the bulge). 

As elliptical galaxies, bulges are diverse and heteroge- 
neous objects. Big bulges of lenticulars and early-type spi- 
rals are similar to low-luminosity elliptical galaxies. On the 
contrary, small bulges of late- type spirals a re reminiscent of 
disks (see the reviews by Korm endvlll993t IWvse etal.|[l997t 
Kormend v & Kenni cutt 2004). Some of them h ave a quite 
complex structure a nd host nuclear rings (see iButal j 1995: 
Com eron et al.l 2010, for a compilation), inne r bars (see | Erwi 
2004 for a list), and embedded disks (e.g., IScorza & Bende 



1995 



Ivan den Bosch et al.lll998l:lPizzella et al]l2002l) . Although 



the kinematical properties of many bulges are well described by 
dynamical models of oblate ellipsoid s which are flattened by 
rotation with little or no anisotropy (Kormendv & Illingworth 



rotation witn little or no anisotropy ( Jvormenav &l illingwortr 
1982; lDavies & IlUngworthlj 1983b lFillmord|1986t ICorsini et al 



1999: lPignatelli et al.ll200U IC appellari et a i1l2006l). the tw isting 
of the bulge isophotes (Lindb ladll956HZaritskv & Lolll986|) and 
the misalignment be t ween the major axes of the bulge and disk 
(IBertola et all I199U IVarela et al.l 119961: iMendez-Abreu et alj 
2008) observed in several galaxies are not possible if the bulge 
and disk are both axisymmetric. These features were interpreted 
as the signature of bulge triaxiality. This ide a is supported by 
the presence of non-circul a r gas motions ( e.g., Gerhard & Vietri 



1986: iBertola et all 1 19891 iGerhard et ail 119891: lBermanll200U 
Falcon-Barroso et alj 20061 : iPizzella et al J 12008]) and a velocity 
gradient along the galaxy minor axis ("e.g.. ICorsini et al] r2003: 
ICoccato et al 12004020051) . 

Perfect axisymmetry is also ruled out when the intrinsic 
shape of bulges is determin ed by statisti c al ana lyses based on 
their observed ellipticities. IBertola et al.l (1 1 99 lb measured the 
bulge ellipticity and the misalignment between the major axes 
of the bulge and disk in 32 SO-Sb galaxies. They found that 
these bulges are triaxial with mean axial ratios (B/A) = 0.86 
and {CIA) = 0.65. (B/A) = 0.79 for the bulges of 35 early- 
type disk galaxies an d (B/A) = 0.71 for the bulges of 35 late- 
type spirals studied by Fathi & Peletierl (120031) . They derived the 
equatorial ellipticity by analyzing the deprojected ellipticity of 
the ellipses fitting the galaxy isophotes within the bulge radius. 
None of the 21 disk g alaxies with morphological types bet ween 
SO and Sab studied bv lNoordermeer & van derHulst(2007) har- 
bors a truly spherical bulge. A mean flattening (C/A) = 0.55 was 
obtained under the assumption of bulge oblateness by compar- 
ing the isophotal ellipticity in the bulge- dominated region to that 
measured in the disk-dominated region. iMosenkov et al.l {2010) 
obtained a median value of the flattening (C/A) = 0.63 for a 
sample of both early and late-type edge-on galaxies in the near 
infrared. They found also that bulges with Sersic index n < 2 can 
be described by triaxial, nearly prolate bulges that are seen from 
different projections, while n > 2 bulges are better represented 
by oblate spheroids with moderate flattening. 

In iMendez- Abreu et ail d2008l hereafter Paper I) we mea- 
sured the structural parameters of a magnitude-limited sample 
of 148 unbarred early-to-intermediate spiral galaxies by a de- 
tailed photometric decomposition of their near-infrared surface- 
brightness distribution. The probability distribution function 
(PDF) of the bulge equatorial ellipticity was derived from the 
distributions of observed ellipticities of bulges and misalign- 
ments between bulges and disks. We proved that about 80% of 
the sample bulges are not oblate but triaxial ellipsoids with a 
mean axial ratio (B/A) - 0.85. The PDF is not significantly de- 
pendent on morphology, light concentration, or luminosity and 



is independent on the possible presence of nuclear bars. This is 
by far the largest sample of bulges studied with this aim. 

In this paper, we introduce a new method to derive the in- 
trinsic shape of bulges under the assumption of triaxiality. This 
statistical analysis is based upon the analytical relations between 
the observed and intrinsic shapes of bulges and their surrounding 
disks and it is applied to the galaxy sample described in Paper 

1. The method was conceived to be completely independent of 
the studied class of objects, and it can be applied whenever tri- 
axial ellipsoids embedded in (or embedding) an axisymmetric 
component are considered. 

The structure of the paper is as follows. The basic description 
of the geometry of the problem and main definitions are given in 
Sect. [2] The statistical analysis of the equatorial ellipticity and 
intrinsic flattening of bulges is presented in Sect. [3] and |4] re- 
spectively. The intrinsic shape of bulges is discussed in Sect. [5] 
The conclusions are presented in Sect. [6] 

2. Basic geometrical considerations 

As in Paper I, we assume that the bulge is a triaxial ellipsoid and 
the disk is circular and lies in the equatorial plane of the bulge. 
Bulge and disk share the same center and polar axis. Therefore, 
the inclination of the polar axis (i.e., the galaxy inclination) and 
the position angle of the line of nodes (i.e., the position angle 
of the galaxy major axis) are directly derived from the observed 
ellipticity and orientation of the disk, respectively. 

We already introduced in Paper I the basic geometrical def- 
initions about the triaxial ellipsoidal bulge and its deprojection 
as a function of the main parameters which describe the prob- 
lem, i.e., the ellipticity e of the projected ellipse, twist angle d 
between its major axis and the line of nodes, galaxy inclination 
9, and orientation <p of the equatorial axes of the bulge with re- 
spect to the line of nodes. However, for the sake of clarity we 
will review again these concepts in this section together with the 
new definitions needed to perform our statistical approach. 

2.1. Direct problem: from ellipsoids to ellipses 

Let (x,y, z) be the Cartesian coordinates with the origin in the 
galaxy center, the jc-axis and y-axis corresponding to the prin- 
cipal equatorial axes of the bulge, and the z-axis corresponding 
to the polar axis. As the equatorial plane of the bulge coincides 
with the equatorial plane of the disk, the z-axis is also the polar 
axis of the disk. If A, B, and C are the lengths of the ellipsoid 
semi-axes, the corresponding equation of the bulge in its own 
reference system is given by 



2 2 2 

xr y z £_ _ 



(1) 



It is worth noting that we do not assume that A > B > C, as 
usually done in the literature. 

Let (x',y',z') be now the Cartesian coordinates of the ob- 
server system. It has its origin in the galaxy center, the polar 
z'-axis along the line of sight (LOS) and pointing toward the 
galaxy. The plane of the sky lies on the (x',y') plane. 

The projection of the disk onto the sky plane is an ellipse 
whose major axis is the line of nodes (LON), i.e., the inter- 
section between the galactic and sky planes. The angle be- 
tween the z-axis and z'-axis corresponds to the inclination of 
the galaxy and therefore of the bulge ellipsoid; it can be de- 
rived as 6 — arccos(ii/c) from the length c and d of the two 
semi-axes of the projected ellipse of the disk. We defined <p 
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(0 < <p < n/2) as the angle between the x-axis and the LON 
on the equatorial plane of the bulge (x,y). Finally, we also de- 
fined tff (0 < iff < n/2) as the angle between the x'-axis and the 
LON on the sky plane (x',y'). The three angles 9, <p, and if/ are 
the usual Euler angles and relate the reference system (x,y,z) 
of the ellipsoid with that (x',y',z') of the observer by means of 
three rotations (see Fig.[TJ. Indeed, because of the location of the 
LON is known, we can choose the x'-axis along it, and conse- 
quently it holds that if/ = 0. By applying these two rotations to 
Eq.Q]it is possible to derive the equation of the ellipsoidal bulge 
in the reference system of the observer, as well as the equation 
of the ellipse correspon ding to its projection on the sky plane 
(Sim onneau et alJ [l998). Now, if we identify the latter with the 
ellipse projected by the observed ellipsoidal bulge, we can de- 
termine the position of its axes of symmetry x e and y e and the 
lengths a and b of the corresponding semi-axes. The x e - axis 
forms an angle 6 with the LON corresponding to the x'-axis of 
the sky plane. We always choose < 6 < n/2 such that a can be 
either the major or the minor semi-axis. If a corresponds to the 
major semi-axis then b is the length of the minor semi-axis. If 
a corresponds to the minor semi-axis then b is the length of the 
major semi-axis. Later in this paper, when we will present our 
statistical analysis we will find that this riddle is solved because 
the two possibilities coincide, and one is the mirror image of the 
other. 

From the previous considerations (see ISimonneau et alj 

Il998l for details) we have that the equations relating the length 
of the semi-axes of the projected ellipse with the length of the 
semi-axes of the intrinsic ellipsoid are given by 



crb 1 = A 2 C 2 sin 2 0cos 2 + fl 2 C 2 shrflsin 2 0+A 2 ,B 2 cos 2 e, 

A 2 (cos 2 <p + cos 2 flsin 2 tp) + B 2 (sin 2 <j> + cos 2 9 cos 2 <p) + C 2 sin 2 
(B 1 -A 2 ) cos 9 sin 20 



'/ + b 1 = 
tan 2a = 



A 2 (cos 2 9 sin 2 <p - cos 2 (/>) + fi 2 (cos 2 9 cos 2 <p - sin 2 <f>) + C 2 sin 2 9 



(2) 
(3) 

(4) 



If the ellipsoidal bulge is not circular in the equatorial plane 
(A + B) then it is possible to observe a twist (6 + 0; see Eq. |4]l 
between the axes of the projected ellipses of the bulge and disk. 

2.2. Inverse problem: from ellipses to ellipsoids 

We will focus now our attention on t he inverse problem, i.e ., 
the problem of deprojection. Followin dSimonneau et alj (|1998). 
from Eqs. [2] [3] and|4] we are able to express the length of the 
bulge semi-axes (A, B, and C) as a function of the length of the 
semi-axes of the projected ellipse (a, b) and the twist angle (6). 

For the sake of clarity, we rewrite here the corresponding 
equations but in a different way with respect to Paper I. First, we 
define 



K 1 



where 



b 2 



[1 + e cos 26] . 



1 < e < 1, 



(5) 



(6) 



is, in some sense, a measure of the ellipticity of the observed 
ellipse. Therefore, K 2 is a positive measurable quantity. 
From Eqs. [2EI and|4]it results 



K 2 = 



A 2 + B 2 



[l+E cos 20] . 



(7) 



where 

A 2 - B 2 



A 2 + B 2 



!<£<!, 



(8) 



measures the intrinsic equatorial ellipticity of the bulge. 

With this notation we can rewrite the equations for the semi- 
axes of the bulge in the form 



A 2 
B 2 
C 2 



K 2 1 + 



K 1 



e sin 28 tan<; 
1 + e cos 25 cos ( 

e sin 26 coti; 
1 + e cos 25 cos ( 
2 e cos 26 



2 e cos 9 sin 28 



sin" 6 ( 1 + e cos 26) sin" 6 (1 + e cos 26) 



(9) 
(10) 

cot 2 0).(11) 



The values of a,b,6, and 6 can be directly obtained from 
observations. Unfortunately, the relation between the intrinsic 
and projected variables also depends on the spatial position of 
the bulge (i.e., on the <p angle), which is actually the unique un- 
known of our problem. Indeed, this will constitute the basis of 
our statistical analysis. 

2.3. Characteristic angles 

There are physical constraints which limit the possible values of 
4>, su ch as the positive length of the three semi-axes of the ellip- 
soid (ISimonneau et al. 1998). Therefore, we define some char- 
acteristic angles which constrain the range of (p. Two different 
possibilities must be taken into account for any value of the ob- 
served variables a, b, 6 and 0. 

The first case corresponds to a > b. It implies that e > from 
Eq.|6]and A > B from Eqs.|9]and[T0] For any value of <p, A 2 > K 2 
and K 2 is always positive according to Eq.|7] On the other hand, 
B 2 and C 2 can be either positive or negative depending on the 
value of <p according to Eqs.fTOlandfTTl respectively. This limits 

d>B- The 



the range of the values of ( 
angle <ps is defined by B 2 - 



. B is positive only for . 
0inEq.[l0]as 



tan0£ 



e sin 26 



(12) 



cos 6 (1 + e cos2<5) 
Likewise, C 2 is positive only for values of (p < (pc- The angle 



be is defined by C 2 = in Eq.fTTIas 

2 e sin 26 cos 9 



tan 20c 



e cos 26 (1 + cos 2 6) - sin" 9 



(13) 



Thus, if a > b then the values of <p can only be in the range 
4>b < </> < <Pc- 

The second case corresponds to a < b. It implies that e < 
(Eq.© and A < B (Eqs.|9]and[T0]l. For any value of </>, B 2 > K 2 
and K 2 is always positive according to Eq. [7] But, A 2 and C 2 
can be either positive or negative depending on the value of <p 
according to Eqs.|9]and[TT] respectively. This limits the range of 
the values of (p. A 2 is positive only for <p < </>a- The angle (f>A is 
defined by A 2 — in Eq.|9]as 



tanc 



cos 9(1 + e cos 26) 
e sin 26 



(14) 



Likewise, C 2 is positive only for values of (p > (pc- The angle 
<pc is given in Eq. [13] Thus, if a < b, then the values <p can only 
be in the range <pc < <p < <Pa- 
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Fig. 1. Schematic three-dimensional view of the ellipsoid geometry. The bulge ellipsoid, the disk plane, and the sky plane are shown 
in red, blue, and green, respectively. The reference systems of the ellipsoid and of the observer as well as the LON are plotted with 
thin solid lines, thin dashed lines, and a thick solid line, respectively. The bulge ellipsoid is shown as seen from an arbitrary viewing 
angle (left panel), along the LOS (central panel), and along the polar axis (i.e., the z-axis; right panel). 



However, the problem is symmetric: the second case, when 
the first semi-axis of the observed ellipse (which is measured 
clockwise from the LON) corresponds to the minor axis (i.e., 
a < b), is the mirror situation of the first case, when the first 
measured semi-axis of the observed ellipse corresponds to the 
major axis (i.e., a > b). In the second case, if we assume the 
angle ji/2 - 5 to define the position of the major semi-axis a of 
the observed ellipse with respect to the LON in the sky plane, 
and it /2 - to define the position of the major semi-axis A of 
the equatorial ellipse of the bulge with respect to the LON in the 
bulge equatorial plane, then we can always consider a > b and 
A > B. Therefore, e > and E > always. This means that 
we have the same mathematical description in both cases: the 
possible values of (p are <pB < (p < (pc with <Pb and (pc defined by 
Eqs. [12] and [13] respectively. Furthermore, we can rewrite Eqs. 
EM andQI] 



From Eq.[16]we obtain that the semi-axis length B is zero for 
<p — 4>b and it increases when (p goes from <pB to (pc- The semi- 
axis length C is zero for (p = (pc and decreases when <p goes from 
4>b to 4>c- There is an intermediate value <Pbc f° r which B 2 = C 2 . 
This angle is given by 



tan (pBc 



tan 5 

cos 9 



(18) 



the corresponding 



For < 

than A 2 . This implies that in this range of < 
triaxial ellipsoid is completely oblate. 

On the other hand, B 2 < A 2 for all possible values of (p. This 
is not the case for C 2 , because it increases when (p decreases. 
Thus, we can define a new angle cp^c for which C 2 - A 2 . This 
angle is given by 



K 2 [1 + tan 4> B tan0] , 



B 2 = K 2 



tan; 



tan< 



i -> cos* 8 r 

C 2 = K-2tm<p B — — [cot 2 (f> - cot 2 (f> c ] . 
sin 6 



(15) 
(16) 

(17) 



where ips and 4>c are given as a function of the observed vari- 
ables a, b, 6, and 8, i.e., they are known functions for each ob- 
served bulge. 

We can always consider A > B as explained before. But, 
we are not imposing any constraint on the length C of the polar 
semi-axis. According to this definition oblate and prolate triaxial 
ellipsoids do not necessarily have an axisymmetric shape. We 
will define a triaxial ellipsoid as completely oblate if C is smaller 
than both A and B (i.e., the polar axis is the shortest axis of the 
ellipsoid). We define a triaxial ellipsoid as completely prolate if 
C is greater than both A and B (i.e., the polar axis is the longest 
axis of the ellipsoid). If the polar axis is the intermediate axis 
we have either a partially oblate or a partially prolate triaxial 
ellipsoid. A further detailed description of all these cases will be 
given at the end of this section. 



tan 4>ac — cos & tan 6. 



(19) 



For (p < <Pac, C 2 > A 2 > B 2 . Therefore, the corresponding 
triaxial ellipsoid is completely prolate. It is important to notice 
here that this case is physically possible only when <Pac > <Pb, 
and the values of <p are within the range of possible values 
(pB ^ <P ^ <Pc- Therefore, we conclude that for any observed 
bulge (i.e., for any set of measured values of a, b, 6, and 6) 
the corresponding triaxial ellipsoid could be always completely 
oblate, while we are not sure that it could be prolate. 

We define the quadratic mean radius of the equatorial ellipse 
of the bulge in order to extensively discuss all the different pos- 
sibilities. 



A 2 + B 2 



K 2 tancps [cot (pB - cot 2 <p] . 



(20) 



which depends only on the unknown position (p. 

Since A 2 > B 2 , A 2 > R 2 > B 2 but there is always a value (p RC 
corresponding to the case C 2 = R 2 



tan 2<Prc = tan 26 ■ 



1 + cos 2 1 
2 cos 6 



(21) 
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The mean equatorial radius allows us to distinguish oblate 
(C 2 < R 2 ) and prolate (C 2 > R 2 ) triaxial ellipsoids. 
Unfortunately, the situation is more complicated and there are 
four different possibilities for the intrinsic shape of the bulge 
ellipsoid. They are sketched in Fig. [2] and can be described as 
follows 

- If <Pac < <Prc < <Pb the triaxial ellipsoid is always oblate (Fig. 
[2] left panel). It is either completely oblate (i.e., A > B > C) 
iiR> B> C (4>bc < 4>< <t>c) or partially oblate if R > C > B 

(<f>B < < <Pbc)- 

- If <Pac < 4>b < 4>rc the triaxial ellipsoid can be either 
oblate or prolate (Fig. [2] central panel). It is either com- 
pletely oblate if R > B > C {ipse < < (pc), or partially 
oblate if R > C > B (0rc < <P < 4>bc), or partially prolate if 
C>R> (0 B < < 4> RC ). 

- If </>b < (f>Ac < <Prc f° ur different possibilities are allowed for 
the triaxial shape of the bulge ellipsoid (Fig. [2] right panel). 
It is either completely oblate if R > B > C (<Pbc < <P < <PcX 
or partially oblate if R > C > B (0#c < < <Pbc), or partially 
prolate if A > C > R (<Pac < <t> < <Pbc), or completely prolate 
(i.e., C> A> B) if C>A>R(<p B <4>< 0ac)- 



3. Equatorial ellipticity of bulges 

In Paper I we focused on the equatorial ellipticity defined in Eq. 
[8] This is a straightforward definition resulting from the equa- 
tions involved in projecting and deprojecting triaxial ellipsoids. 
It allows us to solve the problem of inverting an integral equa- 
tion in order to derive the PDF of the equatorial ellipticity of 
bulges. However, the usual axial ratio B/A is a more intuitive 
choice to describe the equatorial ellipticity of the bulge when 
only one galaxy is considered. Therefore, we redefine the equa- 
torial ellipticity as Z = B 2 /A 2 . Adopting a squared quantity gives 
us the chance of successfully performing an analytic study of the 
problem. By taking into account Eqs.[l5]and[T6] we obtain 



2 sin< 



>B 



z= B 2 = tan (0- B ) = 1 

A 2 tan0 sin0B + sin (20 - <ps) 



Z r = 



Z — for <p = (ps, while the limiting value of Z for 
tan (0 C - </>b) 2 sin B 



tanc 



= 1 - 



sin 4> B + sin (2 <p c - <f> B ) ' 



(22) 

4>c is 
(23) 



When (p is between (pB and <pc, the value of Z reaches a max- 
imum given by 



Zm — 



1 - sin ( 
1 + sin i 



which is observed when (p corresponds to 



n (p B 



(24) 



(25) 



where cpM is always larger than <ps- The value Z decreases for 
4> > <Pm, after reaching its maximum Zm at <p = 0m- Z = for 
<p = n/2. But, it is not necessary to study the behaviour of Z for 
<pc < (p < tt/2 since this range of <p is not physically possible. 

Therefore, as soon as cp increases from (ps to (pc there are 
two possible cases for 0m and the corresponding trend of Z. If 
0c > 0m> the value of Z reaches the maximum Zm for = 0m- 
For larger values of <p it decreases, reaching the limit value Zq for 
= (pc- If 0c < 0m, Z does not reach the maximum value given 



by Eq.[24] In this case, the maximum value of Z corresponds to 

z c . 

We also derive for each observed bulge the mean value (Z) 
of its equatorial ellipticity. From Eq.[22l 



•tpc 



(Z) = 



Z(0) d<p = 



— r 



tan (p B 



b c - 4> B cos (0c - (p B ) sin (p B 



In- 



sin 0c 



(26) 



To perform a more exhaustive statistical analysis, we com- 
pute for each observed bulge the probability P(Z) corresponding 
to < Z < Zc by taking into account that can take any value 



in the range 0# 



< d> < 



be with the same probability given by 



P(0) = 



1 



(27) 



P(Z) = Yj P(4>) \d<p I dZ\, where the sum is defined over all 
the values which solve Eq. [22] The probability P(Z) allows 
us to compute some characteristic values of Z, such as the me- 
dian value Zi/2. It is defined in such a way that the integrated 
probability between Z = and Z\/2 is equal to the integrated 
probability between Z1/2 and Zc- 

The distribution of the sample bulges as a function of their 
maximum, mean, and median equatorial ellipticity is plotted in 
Fig.0 

Moreover, we define the confidence interval (Zi/6,Zs/6) 
where the integrated probability is 67%. The integrated proba- 
bilities between Z = and Z1/6 and between Z = and Z5/6 are 
1/6 and 5/6, respectively. To this aim, we introduce three charac- 
teristic values of in the range between <ps and 0c- According 
to the probability P(<p) given in Eq. [27] they are 



\/2 = ^ + 2^ B ' 



01/6 = + ^<PB, 

5 1 
05/6 = t4>C + t4>B- 



1 



(28) 
(29) 
(30) 



We have seen that Z has a different behaviour for 0c < 0m 
and for 0c > 0m- Therefore, we will separately study these two 
cases in order to derive P(Z) and the corresponding distribution 
of equatorial ellipticities. 

3. 1 . Bulges with <p c < M 

If 0c < 0m, the value ofZmonotonically increases from Z(0b) = 
to Zc = Z(0c). There is only one value of corresponding to 
any given value of Z. Thus the integrated probability P(Z) from 
Z = to Z = Z1/6, Z1/2, and Z5/6 is equal to the integration 
of P((p) from = (ps to = 0^ 6 , 0°^ 2 , and 0^ 6 , respectively. 
Consequently, the median value is 



Z 1/2 = Z(0? /2 ) = 1 



2 sin (ps 
sin <pB + sin 0c' 



and the limits of the confidence interval are 



-1/6 



Z(0i/ 6 ) = 1 



2 sin< 



sin0 B + sin(|0 c + \<pe 



r 



(31) 



(32) 




Fig. 2. The lengths A, B, and C of the semi-axes of the bulge ellipsoid and its mean equatorial radius R as a function of the angle 
(p. The solid lines correspond to the ranges of physically possible values of A, B, C, and R, while the dotted lines show their overall 
trends within < < n/2. A triaxial bulge with 4>ac < </>rc < </>b, 4>ac < 4>b < <Prc, an d <Pb < 4>ac < <Prc is shown in the left, 
central, and right panel, respectively. 



and 



Z 5/6 = Z(^ /6 ) = 1 - 



2 sin<pB 



sin <p B + sin (f 0c - §0 B ) 



P(Z) 



In this case, the probability P(Z) is 
1 sin^B 



h (1 -Z) A /(l-Z) 2 -sin 2 B (l+Z) 2 



which increases monotonically between 
1 



P(0) 
and 



tan <pB, 



P(Z C ) = - 

<Pc- 4>B 4 



> B )] 2 



1 [sin <pB + sin (2(pc - 
sin B cos (2 4> c - <f> B ) 



(33) 



(34) 



(35) 



(36) 



The probability P(Z) given in Eq. [34] strongly peaks at Z = 
Zc in such a way that Z1/2 is close to Zq. For this reason, al- 
though the right portion (Z\/2,Zs/6) of the confidence interval 
(Z1/6, Zg/6) is not large, the confidence interval spans a large frac- 
tion of the total range between and Zc- This is the case for the 



bulge of MCG -02-33-017 (Fig. H top panel). Using the mean 
(Z) and median Z1/2 values to describe the equatorial ellipticity 
of these kinds of bulges is a poor approximation. 



3.2. Bulges with 4> c 

For <pc > 4>m, Z monotonically increases from Z((ps) = to 
Zm = Z(0m) and then it monotonically decreases from Zm to 
Zc = Z((pc). For < Z < Zc there is only one value of <p for 
each value of Z, while for Zc < Z < Zm there are two values of (p 
which correspond to each value of Z. There is a discontinuity in 
P(Z) for Z = Zc, which corresponds to the value <p' c — | - (<pc - 
<Pb)- dZ/d(p = for <p — <Pu and the probability P(Z) becomes 
infinity at Z = Zm- It is not possible to compute directly the 
median value Z1/2 and confidence interval (Zi/^Zs/e) from P((f>) 
in Eq.[27] Therefore, we need to rewrite P(Z) as 



sin 



— \ 2 sin^B 



(1-z) Vd-z) 2 -sin 2 B (l+z) 2 



< Z < Z c , 
Zc < Z < Zm. 



(37) 
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Fig. 3. The distribution of the 148 sample bulges as a function 
of their maximum (top panel), mean (middle panel), and median 
(bottom panel) equatorial ellipticities plotted with a solid line. 
In each panel, the dotted line refers to the distribution of the 115 
sample bulges with <pc > <Pm- 



There are different values for Z1/6, Z1/2 
on whether </>y 2 is smaller or greater than 
to the discontinuity in P(Z). 



and Z5/6 depending 
•' c which corresponds 



Fig. 4. PDF of the equatorial ellipticity for three sample bulges. 
MCG -02-33-017 (top panel) hosts a bulge with <p c < <pu- 
NGC 1 107 (middle panel) hosts a bulge with <pc > (pM and <p^ 2 < 
<p' c . NGC 4789 (bottom panel) hosts a bulge with <pc > <Pm and 
0?/2 > 0c ■ ^ n eacn P ane l> tne vertical line shows the median Z1/2 
value, the arrow corresponds to the maximum value of Z, and 
the hatched area marks the confidence interval (Zi/g, Z5/6) corre- 
sponding to 67% probability. 



For </> 



'„ the values of Z1/2 and Z1/6 are given by Eqs. 
[3T1 and [32] respectively. But, there are two possible values for 



1/2 



Z5/6 depending on the value of 
given by the Eq. [33] If cf>° 5/6 > , 



05/6- If 05/6 < 0C Z 5/6 is 

b' the corresponding values of 
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Z are on the right side of the discontinuity (i.e., two values of Z 
correspond to a given value of (p > (p' c ). In this case 



^5/6 



1 



2 sin( 



sin^B + cos — ( f s 



(38) 



which corresponds to Z(05/6) with (ps/6 = tt/4 — 0c/12 + 70b/12. 



For <p"i 2 > 4>' c the value of Z\ /2 is given by 



2 sin< 



'1/2 



sin<^>B + cos S£Z£° 



(39) 



and it corresponds to Z(<pi/2) with 0i/2 = tt/2 - 0c/2 + 3 0b/2. 
Likewise, Z5/6 is given by Eq. [38] But, for Zi/g we have two 
possibilities according to the value of (p®^. If (p°^ 6 < <p' c then 

Z1/6 is given by Eq.[32] If <p^ 6 > <p' c the corresponding values of 
Z are on the right side of the discontinuity, and it is 



-1/6 



1 



2 sin<pB 



sm<pB + cos 



5 (<pc- 



(40) 



'1/6 



tt/4 



:/12 



which corresponds to Z((p\/(,) with 
ll^j/12. 

For (f>c > <Pm the probability P(Z) in Eq. [37]peaks strongly 
at Zm and therefore the median Z1/2 and maximum Zm values 
of the equatorial ellipticity are very close and the confidence in- 
terval (Zi/6,Z5/6) is narrow. This is the case for the bulges of 
NGC 1 107 (Fig.S middle panel) and NGC 4789 (Fig.0 bottom 
panel). We conclude that for these types of bulges the statistics 
we have presented here are representative of their intrinsic equa- 
torial ellipticity. 



3.3. Statistics of the equatorial ellipticity of bulges 

The distribution of the maximum equatorial ellipticity (corre- 
sponding to either Zq for bulges with (pc < (pM or Zm for bulges 
with (pc > <Pm) peaks at Zm > 0.9 (Fig. [3] top panel). These 
are nearly circular bulges (B/A = 0.95). But, we conclude that 
a large fraction of the sample bulges are strong candidates to be 
triaxial because 41% of them have Z M < 0.80 (B/A < 0.89). 
This result is in agreement with our previous finding in Paper I 
and with the analysis of the distribution of mean (Fig. [3] mid- 
dle panel) and median (Fig. [3] bottom panel) ellipticities. In fact, 
we find that 64% and 53% of our bulges have (Z) < 0.8 and 
Z1/2 < 0.8, respectively. The mean values of (Z) and Z1/2 are 
0.68 and 0.73, respectively. 

The width of the confidence interval (Zi/6,Z5/g) correspond- 
ing to a 67% probability is related to the accuracy of the Z 
measurement. The narrowest confidence intervals are found for 
bulges with (pc — > tt/2 and (ps — > 0. This implies that (p' c — > (p B 
and Zq > Zm- For these bulges the discontinuity in P(Z) is al- 
most negligible. The case with (pc = n/2 and (ps = corre- 
sponds either to spherical bulges (i.e., e — 0) or to bulges with 
a circular equatorial section (i.e., tan 25 = 0). Consequently, the 
bulges with B « A are among those characterized by the nar- 
rower confidence interval and better determination of Z. We can 
select all sample objects for which the Z measurement is only 
slightly uncertain. They are the 115 galaxies with (pc > <Pm- 
The distribution of these selected bulges as a function of their 
Zm, (Z), and Z1/2 is plotted in Fig.[3]too. The fraction of bulges 
with Zm < 0.8 is 33%. It is significantly smaller than the 41% 
found for the complete sample, because the selected sample is 
biased toward bulges with B w A including all the bulges with 



a circular (or nearly circular) equatorial section. The fraction of 
selected bulges with <Z> < 0.8 and Z 1/2 < 0.8 is 55% and 43%, 
respectively. 

4. Intrinsic flattening of bulges 

The axial ratio C/A usually describes the intrinsic flattening F 
of a triaxial ellipsoid if A > B > C. Since we have no constraints 
on the lengths A, B, and C, we redefine the flattening as 



, s C 2 2C 2 

FW = F = a^7W> 



(41) 



by using the lengths C and R of the polar semi-axis and the mean 
equatorial radius given by Eqs.[TTIandl2"0l respectively. 



F(<p) = F e 



cot 2<p - cot 2 
cot<pB - cot I 



sin <pB sin {2<pc - 2<p) 
sin 2<pc sin (2<p - <ps) 



where 



F = 



2 cos 2 6 
sin 2 6> 



(42) 



(43) 



accounts for the effect of inclination. The angle also enters in 
the definition of the two angles <ps and <pc in Eqs.fT2land[T3l re- 
spectively. Adopting a squared quantity for F gives us the chance 
of successfully performing an analytic study of the problem as 
was done for the equatorial ellipticity Z in Eq. [22] 

Since dF(<p)/d(p < 0, the function F(<p) is monotonically de- 
creasing with a maximum Fm &t(f> = <pB given by 



F m = F, 



sin (2 4> c - 2 <p B ) 
sin 2 d>c 



(44) 



If increases from $b to (pc, the value of F((p) decreases to 
zero at <p — (pc- According to Eq.|44] for Fm < 1 the triaxial el- 
lipsoids are oblate, with some of them being partially oblate and 
others completely oblate. For Fm > 1 the triaxial ellipsoids can 
also be partially prolate and in some extreme cases completely 
prolate. 

From Eq.[42] we compute the mean value (F) of the intrinsic 
flattening as 



(F) = 



(p c -, 



J<f>„ 



= Ft 



sin <ps 
' sin 2 (pc 
■ cos (2 (pc 



F((p) d(p = 

sin (2 (p c - (p B ) , I sin (2 (p c - (p B ) ) 
2 ((p c - (p B ) \ sin (p B j 
- <Pb)] ■ 



(45) 



Since F(<p) is a monotonic function (i.e., each value of F 
corresponds to only one value of (p), the integrated probability 
P(F) between F(<pc) = and some characteristic value F* = 
F(0») is equal to the integral of P((p) between 0» and (pc- Then, 
it is straightforward to compute the median value F1/2 of the 
intrinsic flattening which corresponds to the median value <p^ 2 = 
((Pc + (f> B )/2 



F1/2 = Ffl 



sin (p B sin ((p c - <f> B ) 
sin (2 (p c ) sin (p c 



(46) 



As was done for the equatorial ellipticity, we can define also 
for the flattening a confidence interval (Fj/e, F5/6) where the 
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integrated probability is 67%. In fact, the integrated probabilities 
between F — and F\/§ and between F — and F5/6 are 1/6 and 
5/6, respectively. We have 



'1/6 



sin <pB 



sin 



sin (2 C ) sin( 



(47) 



which corresponds to ( ^ 6 given in Eq.[29j and 

_ _ r sin0B sin(§0c-§0s) 

^5/6 - re — — T, n — \' 

sin (2 0c) sin (l0c + |0 B ) 



(48) 



which corresponds to 0j? , g given in Eq. [30] The distribution of 
the sample bulges as a function of their maximum, mean, and 
median intrinsic flattening is plotted in Fig. [5] 

It is possible to perform a more exhaustive statistical analysis 
by defining the probability P(F) of having a flattening F as 



P(F) = k 



1 



A F 2 +B F + Co 



where 



ko = 



cos 2 6 sin (20c - <f>B) 
sin 2 8 (<pc - <Pb) sin 20c sin 4>b 
1 



sin (pB 

4 cos 2 6 cos (20c - 0b) 
sin 2 sin 20c sin0# 



C„ = 



4cos 4 6 



sin 4 sin 2 20c 



(49) 

(50) 

(51) 

(52) 

(53) 
(54) 



where ko, Aq, and Co are always positive, while Bo > for 20c - 
B < 7r/2 (0c < M ) and B < for 20 c - 0b > n/2 (0 C > 0m). 
All these quantities can be computed directly for each observed 
bulge, indeed they depend only on the measured values of a, b, 
6, and 9 through the angles <pg and 0c- 

In Sect. l3.ll we found that the confidence interval (Z1/6, Z5/6) 
of equatorial ellipticity for a bulge with 0c < 0m is wide. For 
this reason, the median Z1/2 and mean (Z) values are not repre- 
sentative of the equatorial ellipticity of the bulge. The same is 
true for (F1/6, Fg/6) because the probability function P(F) peaks 
at F — and slowly decreases as soon as F increases. As a 
consequence, the median F1/2 and mean (F) values are not rep- 
resentative of the intrinsic flattening of the bulge. This is the case 
for the bulge of MCG -02-33-017 (Fig.© left panels) 

On the contrary, if 0c > 0m then Bo < 0, and the probability 
function P(F) peaks at the most probable value 



1 B 

Fmp = -2A ' 
and it quickly decreases to 
k) 



P(0) 



Co' 



(55) 



(56) 



and to zero for F < Fmp and F > Fmp, respectively. The con- 
fidence interval (F1/6, F5/6) is narrow. The median F1/2, mean 
(F), and the most probable value Fmp are close to each other 
and all of them are representative of the intrinsic flattening. This 
is the case for the bulge of NGC 4789 (Fig.© right panels). 
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Fig. 5. The distribution of the 148 sample bulges as a function 
of their maximum (top panel), mean (middle panel), and median 
(bottom panel) intrinsic flattening, plotted with a solid line. In 
each panel, the dotted line refers to the distribution of the 115 
sample bulges with 0c > 0m- 



4. 1 . Statistics of the intrinsic flattening of bulges 

The distribution of the maximum intrinsic flattening (Fig.© top 
panel) shows that 12% of the sample bulges have Fm < 1 (i.e., 
they are either completely or partially oblate triaxial ellipsoids). 
Judging by Fm, the majority of sample bulges could be highly 
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MCG -02-33-017 NGC 4789 




Fig. 6. The intrinsic flattening as a function of the angle <p (top panels) and its PDF (bottom panels) for two sample bulges. MCG - 
02-33-017 (left panels) and NGC 4789 (right panels) host a bulge with Bq > and Z?o < 0, respectively. In the bottom panels, the 
vertical line shows the median F1/2 value and the hatched area marks the confidence interval (F1/6, corresponding to 67% of 
probability. 



elongated along the polar axis. However, these highly elongated 
bulges are not common. Indeed, after excluding from the com- 
plete sample the bulges with Fm < 1, only 19% (18% if we con- 
sider only the selected sample of 115 bulges) of the remaining 
bulges have a probability greater than 50% to have an intrinsic 
flattening F > 1 and there are no bulges with more than a 90% 
probability of having F > 1 (Fig. |7J. This is in agreement with 
the results based on the analysis of the distribution of the mean 
(Fig EJ middle panel) and median (Fig. [5] bottom panel) intrinsic 
flattening. We find that 78% of the sample bulges have (F) < 1, 
and 83% have F1/2 < 1. They are oblate triaxial ellipsoids. 

The large number of sample bulges with Fm > 1 with re- 
spect to those which are actually elongated along the polar axis 
is due to a projection effect of the triaxial ellipsoids. For any <p 
the contribution of inclination 6 to the value of F is given by 
Fg as defined in Eq. [43] However, the intrinsic flattening scales 
with Fg, whereas the probability P(F) scales with l/Fg. Thus, 
the probability to have large Fg values (and large Fm values) is 
very small. For instance, the probability to have the maximum 
Fm value given by Eq.l44lis 



P(F M ) = 



1 



1 sin^B sin 2 0c 



2 ((f>c ~ <Pb) Fg sin (2 <p c - 4>b) 



(57) 



We conclude that Fm is not a good proxy for the intrinsic flat- 
tening of a bulge, although Zm is a good proxy for equatorial 
ellipticity. 

The distribution of the selected bulges with <pc > 0m as a 
function of their Fm, (F), and F1/2 is also plotted in Fig. [5] The 
fraction of oblate triaxial ellipsoids is rather similar to that of 
the complete sample, being 10%, 78%, and 83% if we consider 
bulges with Fm < 1, (F) < 1, and F[/2 < 1, respectively. The 
mean values of (F) and F1/2 are 0.88 and 0.71, respectively, for 
the complete sample, and 0.86 and 0.75, respectively, for the 
selected sample. 



5. Intrinsic shape of bulges 

The distributions of the equatorial ellipticity and intrinsic flat- 
tening of bulges have been studied in Sects. [3] and [4] as two in- 
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dependent and not correlated statistics. It is possible to find the 
relation between them from Eqs.lBlandflTI 



E 2 



Y u sin 2<pc + sin <ps cos (20c - <Pb) 



sin <ps 



sin (2 4> c - 4>b) 



(58) 



to constrain the intrinsic shape of an observed bulge with the 
help of the known characteristic angles 4>b and <pc, which depend 
only on the measured values of a, b, 5, and 9. Eq. [58] can be 
rewritten as a function of the axial ratios B/A and C/A as 



2 sin (20 c ) C 2 



A 2 



sin( 




sin q>B 



l + A 2 
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B 2 

sin d> B cos (20c - cf> B ) 1 1 + — 



(59) 



MCG -02-33-017 



Since B/A and C/A are both functions of the same variable 
0, their probabilities are equivalent (i.e., for a given value of 
B/A with probability P(B/A), the corresponding value of C/A 
obtained by Eq. g9]has a probability P(C/A) = P(B/A)). This 
allows us to obtain the range of possible values of B/A and C/A 
for an observed bulge and to constrain its most probable intrin- 
sic shape by adopting the probabilities P(Z) and P(F) derived in 
Sects. [3]andm respectively. 

An example of the application of Eq. [59] to two bulges of 
our sample is shown in Fig. [8] where the hatched area marks the 
confidence region which encloses 67% of the total probability 
for all the possible values of B/A and C/A. The intrinsic shape 
of bulges with 0c < 0m is less constrained, since the median 
values of B/A and C/A are less representative of their actual val- 
ues. This is the case for the bulge of MCG -02-33-017 (Fig. [H 
top panel). On the contrary, the intrinsic shape of bulges with 
<t>c > 0m is better constrained. This is the case for the bulge of 
NGC 4789 (Fig. [8] bottom panel). 

5. 1 . Statistics of the intrinsic shape of bulges 

Following the above prescriptions, we calculated the axial ratios 
B/A and C/A and their confidence intervals for all the sample 
bulges. There is no correlation between B/A and C/A (Pig. |9j, 
unless only bulges with 0c > 0m are taken into account. The 
range of C/A values corresponding to a given B/A decreases 
as B/A ranges from 1.0 to 0.5, giving a triangular shape to the 
distribution of allowed axial ratios. Circular and nearly circular 
bulges can have either an axisymmetric oblate or an axisymmet- 
ric prolate or a spherical shape. More elliptical bulges are more 
elongated along their polar axis. 

We der i ved t he triaxiality parameter, as defined by 
iFranx et all (Il99ll) . for the 115 sample bulges with a well- 
constrained intrinsic shape (i.e., those with <p c > 0m) 



1 



T = 



(I) 2 



MS 



(60) 



where A, B, and C are the lengths of the longest, intermediate, 
and shortest semi-axes of the triaxial ellipsoid, respectively (i.e., 
A > B > C). This notation is different with respect to that we 
adopted in the previous sections. Oblate triaxial (or axisymmet- 
ric) ellipsoids can be flattened either along the y-axis on the 
equatorial plane of the galaxy or along the polar axis. Prolate tri- 
axial (or axisymmetric) ellipsoids can be elongated either along 
the jc-axis on the equatorial plane of the galaxy or along the 
polar axis. Therefore, prolate bulges do either lie on the disk 
plane (and are similar to bars) or do stick out from the disk (and 
are elongated perpendicularly to it). This change of notation is 
needed to compare our results with those available in literature. 

The triaxiality parameter for bulges with 0c > 0m is char- 
acterized by a bimodal distribution (Fig.fTOb with a minimum at 
T = 0.55 and two maxima at T = 0.05 and T - 0.85, respec- 
tively. According to this distribution, 65% ± 4% of the selected 
bulges are oblate triaxial (or axisymmetric) ellipsoids (T < 0.55) 
and the remaining 35% ±4% are prolate (or axisymmetric) triax- 
ial ellipsoids (T > 0.55). The uncertainties for the percentages 
were estimated by means of Montecarlo simulations. Since T is 
a function of 0, we generated 10000 random values of in the 
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Fig. 8. Relation between the axial ratios B/A and C/A for two 
sample bulges. MCG -02-33-017 (upper panel) hosts a bulge 
with <pM < <pc an d NGC 4789 (lower panel) hosts a bulge with 
4>m > <Pc- The probability associated with each value of B/A and 
its corresponding value of C/A (thick solid line), the value of 
C/A as a function of B/A (dotted line), the maximum value of 
the equatorial ellipticity (arrow), the median values of B/A (ver- 
tical thin solid line) and C/A (horizontal thin solid line), and the 
confidence region which encloses all the possible values of B/A 
and C/A within a 67% probability (hatched area) are shown in 
both panels. 



range between 0# and 0c for each bulge and derived the corre- 
sponding distributions of B/A and C/A according to their PDFs. 
From B/A and C/A we calculated the distribution of T and its 
standard deviation, which we adopted as uncertainty. 

We investigated the cause of such a bimodality by separat- 
ing the bulges according to their Sersic index (n) and bulge-to- 
total luminosity ratio (B/T). Both quantities were derived for 
each sample bulge in Paper I. The Sersic index is a shape pa- 
rameter describing the curvature of the surface-brightness pro- 
file of the bulge. A profile with n — 1 corresponds to an ex- 
ponential law, while a profile with n — 4 corresponds to an 
r l l 4 law. The bimodality is driven by bulges with Sersic index 
n > 2 (Fig.fTOl upper panel), or alternatively, by bulges of galax- 
ies with B/T > 0.3 (Fig. [10] lower panel). In fact, the sample 
of bulges with 0c > <Pm and the two subsamples of bulges with 
n > 2 and of bulges in galaxies with B/T > 0.3 are character- 
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ized by the same distribution of T, as confirmed at high con- 
fidence level (> 99%) by a Kolmogorov-Smirnov test. We find 
that 66%±4% of bulges with n > 2 have T < 0.55. Their number 
decreases as T increases from to 0.55. The remaining bulges 
have T > 0.55 and their number increases as T ranges from 

0. 55 to 1. A similar distribution is observed for the bulges of 
galaxies with B/T > 0.3. 67% ± 4% of them host a bulge with 
T < 0.55. Instead, the distribution of the triaxiality parameter of 
bulges of galaxies with B/T < 0.3 is almost constant with a peak 
at T = 0.05. This is true also for the bulges with n < 2, although 
to a lesser degree. 

The two subsamples of bulges with n < 2 and n > 2 are dif- 
ferent, as confirmed by a Kolmogorov-Smirnov test (99% con- 
fidence level). In particular, the fraction of oblate axisymmetric 
(or nearly axisymmetric) bulges (T < 0.1) is remarkably higher 
for n < 2 (27% ± 4%) than for n > 2 (14% + 3%). The fraction 
of triaxial bulges (0.1 <T< 0.9) is lower for n < 2 (7 1 % ± 5%) 
than for n > 2 (76% ± 3%). The fraction of prolate axisymmetric 
(or nearly axisymmetric) bulges (T > 0.9) for n < 2 is 2% + 2%, 
but 11% + 3% for n> 2. 

The two subsamples of bulges of galaxies with B/T > 0.3 
and B/T < 0.3 are different too, as confirmed by a Kolmogorov- 
Smirnov test (99% confidence level). The fraction of oblate ax- 
isymmetric bulges (T < 0.1) is significantly higher for bulges 
of galaxies with B/T < 0.3 (22% + 4%) than for B/T > 0.3 
(16% + 2%). The fraction of triaxial bulges (0.1 < T < 0.9) is 
significantly lower for B/T < 0.3 (67% ±4%) than for B/ T > 0.3 
(78% ± 3%). A few prolate bulges (T > 0.9) are observed for 
B/T < 0.3 (11% ± 3%) and B/T > 0.3 (6% ± 2%). The distribu- 
tion of bulges with n < 2 and bulges of galaxies with B/T < 0.3 
appears to be the same at a high confidence level (> 99%) as 
confirmed by a Kolmogorov-Smirnov test. 

Bulges with (pc > </>m can be divided into two classes: those 
with n < 2 (or B/T < 0.3) and those with n > 2 (or B/T > 0.3). 
About 70% of bulges with n < 2 are hosted by galaxies with 
B/T < 0.3. The same is true for bulges with n > 2 which are 
mostly hosted by galaxies with B/T > 0.3. This agrees with the 
correlation between n and B/T. 

In order to understand whether the intrinsic shape is corre- 
lated with some of the bulge properties we measured in Paper 

1, we plotted the axial ratios C/A and B/A and the triaxiality of 
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B/T > 0.3 (dashed line) are shown in the lower panel. 



the sample bulges with 4>c > <Pm as a function of their Sersic 
index, /-band luminosity, and central velocity dispersion (Fig. 
[TTI i. As we found in Paper I for the intrinsic equatorial elliptic- 
ity, there are no statistically significant correlations between the 
bulge shape and the bulge Sersic index, luminosity or velocity 
dispersion as pointed out by the low Spearman rank correlation 
coefficient (Fig. [TTJ. However, this could be a selection effect 
since the sample of observed bulges spans over a limited range 
of Hubble types (SO-Sb). 

5.2. The influence of nuclear bars on the intrinsic shape of 
bulges 

Our sample galaxies were selected to not host large-scale bars. 
We checked for their presence in Paper I by a visual inspection 
of both the original image and the residual image we obtained 
after subtracting the best-fitting photometric model. However, 
these selection criteria did not account for the presence of unre- 
solved nuclear bars. Nuclear bars are more elongated than their 
host bulges and have random orientations, therefore they could 
affect the measurement of the structural parameters of bulges and 
consequently of their intrinsic shape. 

In Paper I we built up a set of 1000 artificial images with 
a Sersic bulge, an exponential disk, and a Ferrers nuclear bar 
to study the effects of the bar on the measurements of the pho- 
tometric parameters of bulge and disk. The mean errors on the 
fitted axial ratio and position angle of the bulge ({Aq^), (APAb)) 
and disk ((Agj), (APAa)) and their standard deviations (dAg^, 
6 APAb, dAqd, (5APAd) are given in Table 2 of Paper I. 
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In the present paper, we tested whether including a nuclear 
bar affects the T distribution. For each galaxy, we randomly gen- 
erated a series of 1000 values of q\„ q^, PAb, and PAj. To assess 
whether the bulges appear elongated and twisted with respect 
to the disk due to the presence of a nuclear bar, we assumed 
that the axial ratios were normally distributed around the values 
<7b + (A<?b) and q& + (Aqa) with standard deviations 8Aq\, and 
6Aqa, respectively, and that the position angles were normally 
distributed around the values PAb + (APAb) and PAd ± (Aqa) with 
standard deviations £APAb and £APAd, respectively. We chose 
the PA values that gave the smallest 5 with respect to the ob- 
served one. 

If we assume that all the artificial bulges host a nuclear bar, 
we still obtain a bimodal distribution of T (Fig. PTSb . However, 
the fraction of oblate axisymmetric (or nearly axisymmetric) 
bulges (T < 0.1) is higher (23% + 2%) with respect to the 
observed 18% + 3%. For a more real istic fraction of galaxies 
which host a nuclear bar (i.e, 30%, see lLaine et al.ll20 02: Erwin 
12004 . the resulting distribution of T is consistent within errors 
with the distribution derived in Sect. 15. 1 K Fig.lT2T>. We found that 
20% + 2%, 71% + 3%, and 9% ± 2% of the sample bulges are 
oblate axisymmetric (T < 0.1), triaxial (0.1 < T < 0.9), and 
prolate axisymmetric (T > 0.9), respectively, with respect to the 
18% ± 3%, 74% + 4%, and 8% + 2% previously found. In addi- 



tion, we have also tested the effects of not consider a distribution 
of bar parameters but only the stronger bar included in the simu- 
lations (0.8 x r e , qb = 0.2, and Lb ar = 0.02 x L tot ), i.e., the worst- 
case scenario . If we assume that 30% of our galaxies host this 
kind of nuclear bars the results change strongly, obtaining that 
only 56% ±4% of the sample bulges are triaxial (0.1 < T < 0.9) 
with respect to the 74% ±4% previously found. Considering that 
all galaxies host such kind of nuclear bars the fraction of triaxial 
bulges is 30% ± 3%. 

The measured ellipticity and bulge misalignment with the 
disk of the artificial galaxies without a nuclear bar are smaller 
with respect to the actual values measured for the sample bulges. 
This sets an upper limit to the axisymmetry of the bulges. 

6. Conclusions 

In this work, we have developed a new method to derive the 
intrinsic shape of bulges. It is based upon the geometrical re- 
lationships between the observed and intrinsic shapes of bulges 
and their surrounding disks. We assumed that bulges are triax- 
ial ellipsoids with semi-axes of length A and B in the equatorial 
plane and C along the polar axis. The bulge shares the same cen- 
ter and polar axis of its disk, which is circular and lies on the 
equatorial plane of the bulge. The intrinsic shape of the bulge is 
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recovered from photometric data only. They include the lengths 
a and b of the two semi-major axes of the ellipse, corresponding 
to the two-dimensional projection of the bulge, the twist angle 6 
between the bulge major axis and the galaxy line of nodes, and 
the galaxy inclination 6. The method is completely independent 
of the studied class of objects, and it can be applied whenever a 
triaxial ellipsoid embedded in (or embedding) an axisymmetric 
component is considered. 

We analyzed the magnitude-limited sample of 148 unbarred 
SO-Sb galaxies, for which we have derived (Paper I) the struc- 
tural parameters of bulges and disks by a detailed photometric 
decomposition of their near-infrared surface-brightness distribu- 
tion. 

From the study of the equatorial ellipticity Z = B 2 /A 2 , we 
found that there is a combination of the characteristic angles for 
which the intrinsic shape can be more confidently constrained. 
This allowed us to select a qualified subsample of 115 galax- 
ies with a narrow confidence interval (corresponding to 67% of 
probability) of Z. For example, bulges with B as A are among 
those characterized by the narrower confidence interval and the 
best determination of Z. The fraction of selected bulges with a 
maximum equatorial ellipticity Zm < 0.80 (B/A < 0.89), mean 
equatorial ellipticity (Z) < 0.80 and a median equatorial elliptic- 
ity Z1/2 < 0.80 is 33%, 55% and 43%, respectively. We conclude 
that not all the selected bulges have a circular (or nearly circular) 
section, but a significant fraction of them is characterized by an 
elliptical section. These bulges are strong candidates to be tri- 
axial. In spite of the lower fraction of bulges with a maximum 
equatorial ellipticity smaller than 0.8, Zm is a good proxy for the 
equatorial ellipticity because the selected sample contains all the 
bulges with B m A. 



The analysis of the intrinsic flattening F = 2C 2 /(A 2 + B 2 ) 
shows that only a few bulges of the selected sample are pro- 
late triaxial ellipsoids. Only 22% and 17% have a mean intrinsic 
flattening (F) > 1 or a median intrinsic flattening Ft/2 > 1. re- 
spectively. The fraction rises to 90% when a maximum intrinsic 
flattening Fm > 1 is considered. However, this is due to the pro- 
jection effect of triaxial ellipsoids. Indeed, the fraction of bulges 
which are actually elongated along the polar axis is very small: 
only 18% of bulges with Fm > 1 have a probability greater than 
50% to have an intrinsic flattening F > 1, and there are no bulges 
with more than a 90% probability of having F > 1. Thus, Fm is 
not a good proxy for the intrinsic flattening. 

After considering the equatorial ellipticity and intrinsic flat- 
tening as independent parameters, we derived the relation among 
them in order to calculate for each sample bulge both axial ratios, 
B/A and C/A, and their confidence intervals. As already found 
for Z and F, the axial ratios are better constrained for the selected 
sample of 115 bulges. We de rived the triaxiality parameter, as 
defined bv lFranx etafl (Il99ll) . for all of them. We found that it 
follows a bimodal distribution with a minimum at T = 0.55 and 
two maxima at T — 0.05 (corresponding to oblate axisymmet- 
ric or nearly axisymmetric ellipsoids) and T = 0.85 (strongly 
prolate triaxial ellipsoids), respectively. According to this distri- 
bution, 65% of the selected bulges are oblate triaxial (or axisym- 
metric) ellipsoids (T < 0.55) and the remaining 35% are prolate 
triaxial (or axisymmetric) ellipsoids (T > 0.55). This bimodality 
is driven by bulges with Sersic index n > 2 or alternatively by 
bulges of galaxies with a bulge-to-total ratio B/T > 0.3. Bulges 
with n < 2 and bulges of galaxies with B/T < 0.3 follow a simi- 
lar distribution, which is different from that of bulges with n > 2 
and bulges of galaxies with B/T > 0.3. In particular, the sam- 
ple of bulges with n < 2 and the sample of bulges of galaxies 
with B/T < 0.3 show a larger fraction of oblate axisymmetric 
(or nearly axisymmetric) bulges (T < 0.1), a smaller fraction of 
triaxial bulges (0.1 < T < 0.9), and fewer prolate axisymmetric 
(or nearly axisymmetric) bulges (T > 0.9) with respect to the 
corresponding sample of bulges with n > 2 and the sample of 
bulges of galaxies with B/T > 0.3, respectively. 

The different distribution of the intrinsic shapes of bulges 
according to their Sersic index gives further support to the pres- 
ence of two bulge populations with different structural proper- 
ties: the classical bulges, which are characterized by n > 2 and 
are similar to low-luminosity elliptical galaxies, and pseudob- 
ulges, with n < 2 and charact erized by disk-like properties (see 
iKormendv & Kennicut j|2004l for a review). The correlation be- 
tween the intrinsic shape of bulges with n < 2 and those in galax- 
ies with B/T < 0.3 and between bulges with n > 2 and those in 
galaxies with B/T > 0.3 agrees with the correlation between 
the bulge Sersic inde x and bulge-to-total rat io of the host galaxy, 
as rec ently found bv lDrorv & Fished d2007l) and lFisher & Drorvl 
d2008h . 

No statistically significant correlations have been found be- 
tween the intrinsic shape of bulges and bulge luminosity or ve- 
locity dispersion. However, this could be a selection effect since 
the sample bulges span a limited range of Hubble types (SO-Sb). 

The observed bimodal distribution of the triaxiality param- 
eter can be compared to the prop erties predicted b y numerical 
simulations of spheroid formation. ICox et aD d2006l) studied the 
structure of spheroidal remnants formed from major dissipation- 
less and dissipational mergers of disk galaxies. Dissipationless 
remnants are triaxial with a tendency to be more prolate, whereas 
dissipational remnants are triaxial and tend be much closer to 
oblate. This result is consistent with pre vious studies of dissipa- 
tionless and dissipational mergers (e.g.. iBarnesll 1 992t iHern auist 
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ll992tlSnringell2000l:lGonzalez-Garcfa & Balcellsll2005h . In ad- 
alJOOlol) ' used semi-empirical models to pre- 



dition, iHopkins eti 
diet galaxy merger rates and contributions to bulge growth as 
functions of merger mass, redshift, and mass ratio. They found 
that high B/T systems tend to form in major mergers, whereas 
low B/T systems tend to form from minor mergers. In this 
framework, bulges with n < 2, which shows a high fraction of 
oblate axisymmetric (or nearly axisymmetric) shapes and have 
B/T < 0.3, could be the result of dissipational minor mergers. A 
more complex scenario including both major dissipational and 
dissipationless mergers is required to explain the variety of in- 
trinsic shapes found for bulges with n > 2 and B/T > 0.3. 

On the other hand, depending on the initial conditions (see 
IVietri|[l990l and references therein), the final shape of the early 
protogalaxies could also be triaxial. However, high-resolution 
numerical simulations in a cosmologically motivated framework 
that resolves the bulge structure are still lacking. The compari- 
son of a larger sample of bulges with a measured intrinsic shape 
and covering the entire Hubble sequence with these numerical 
experiments is the next logical step in addressing the issue of 
bulge formation. 
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